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Genome-wide Study of Gene Copy Numbers, 
Transcripts, and Protein Levels in Pairs of 
Non-invasive and Invasive Human Transitional 
Cell Carcinomas* 

Torben F. 0rntoft£§, Thpmas ThykjaeiH, Frederic M. Wa!dman||, Hans Wolf**, 
and Julio E. Celis# 



Gain and loss of chromosomal material is characteristic 
of bladder cancer* as well as malignant transformation In 
general The consequences of these changes at both the 
transcription and translation levels Is at present unknown 
partly because of technical limitations. Here we have at- 
tempted to address this question In pairs of non-invasive 
and Invasive human bladder tumors using a combination 
of technology that Included comparative genomic hybrid- 
ization, high density oligonucleotide array-based monitor- 
ing of transcript levels (5600 genes), and high resolution 
two-dimensional gei electrophoresis. The results showed 
that there is a gene dosage effect that in some cases 
superimposes on other regulatory mechanisms. This ef- . 
feet depended (p < 0.015* on the magnitude of the com- 
parative genomic hybridization change, in general (18 of 
23 cases), chromosomal areas with more than 2-foldgain 
of DNA showed a corresponding increase in mRNA tran- 
scripts. Areas with loss of DNA, on the other hand, 
showed either reduced or unaltered transcript levels, Be- 
cause most proteins resolved by two-dimensional gels 
are unknown ft was only possible to compare mRNA and 
protein alterations in relatively few cases of well focused 
abundant proteins, ggtth few exceptions we found a good 
correlation (p < 0.005) between transcript alterations and 
protein levels. The Implications, as wen as limitations, 
of the. approach are discussed. Molecular & Cellular 
Proteomlcs 1&7-45, 2002. 

, Aneupioidy is a common feature of most human cancers 
(1), but littie is known about the genome-wide effect of this 
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phenomenon at both the transcription and translation Jeveis. 
High throughput array studies of the breast cancer cell line 
BT474 has suggested that there Is a correlation between 
DNA copy numbers and gene expression in highly amplified 
areas (2)* and studies of individual genes in solid tumors 
have revealed a good correlation between gene dose and 
mRNA or protein levels in the case of e-erb-82, cyc//r> off , " 
emsf , and N-myc (3-5). However, a high cyclln D1 protein 
expression has been observed, without simultaneous am- 
plification (4) v and a low level of c-myc copy number In- 
crease was observed without concomitant c-myc protein 
overexpresslon (6). . 

(n human bladder tumors, karyotyping, fluorescent In situ 
hybridization, arid comparative genomic hybridization (CGH) 1 
have revealed chromosomal aberrations that seem to be 
characteristic of certain stages of disease progression.. In the 
case of non-Invasive pTa transitional cell carcinomas (TCCs), 
this includes loss of chromosome 9 or parts of it, as well as 
loss of Y in males. In minimally Invasive pTI TCCs, the fol- 
lowing alterations have been reported: 2q-, 11p-, 1q+, 
11q13+. 17q+, ana" 20q+ (7-12). It has bean suggested that 
these regions .harbor tumor suppressor genes and onco- 
genes: however, the large chromosomal areas involved often 
contain many genes, making meaningful predictions of the 
functional consequences of losses and gains yen/ difficult 

In tills Investigation we-have combined gebbme~wide tech- 
nology for detecting genomic gains and losses (CQH) with 
gene expression profiling techniques (mfcroarrays and pro- 
teomlcs) to determine the effect of gene copy number on 
transcript and protein levels In pairs of non-Invasive and In- 
vasive human bladder TCCs, 

EXPERIMENTAL PROCEDURES 

Material— Bladder tumor biopsies, were sampled after informed 
consent was obtained and after removal of tissue for routine pathol- 
ogy examination. By light microscopy tumors 335 and 532 were 
staged by an experienced pathologist as pTa (superficial papillary), 



1 The abbreviations used are: CGH, comparative genomic hybrid- 
ization; TCC, transitional cell carcinoma; LOH, loss of heterozygosity; 
PA-FABP, psoriasis-associated fatty acid-binding protein; 2D, 
two-dimensional. 
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Ba 1 . DNA copy number- and mRNA expression level. Shown from left to r$M are chromosome (Cfir.), GGH profiles, gene location and 
-expression level of specific genes., and overall expression level along the chromosome. A, expression of mRNA In Invasive tumor 733 as 
compared with the non-Invasive counterpart tumor 335. fl, expression of mRNA In Invasive tumor 827 compared with the non-invasive 
^m m SP a ll tUrTl0r 532 ' 1116 averaflS fluorescent signal ratio between tumor DNA and normal DNA Is shown along the ten gth of the chromosome 
-ifeft). The bold curve In the ratio profile represents a mean of four chromosomes and is surrounded by thin curves Indicating one standard 
devtoOoa The central vertfcaf fine $roken) indicates a ratio value of 1 (no change), and the vertical lines next to It (dotted) indicate a ratio of 
0.5 (/eft) and 2.0 fight}, m chromosomes where the non-Invasive tumor 335 used for comparison showed alterations In DNA content, the ratio 
profile of that chromosome is shown to the right of the invasive tumor profile. The colored oars represents one gene each, Identified by the 
running numbers above the bars (the name of the gene can be seen at www.MDLDKfedata.htm|). The bars Indicate the purported location of 
the gene, and the coAvs indicate the expression level of the gene In the invasive tumor compared with the noninvasive counterpart; >2-foid 
Increase plack). >2-fold decrease [blue), no significant change [orange). The bar to the tar right, entitled Expression shows the resulting change 
* tn expression along the chromosome; the colors Indicate that at least half of -the genes wereup-fegulated (Mac*), at least half of the genes 
dowmragulated (blue), or more than half of the genes are unchanged (orange). If a gene was absent in one of the samples and presenTIri" 
another. It was regarded as more than a 2-fold change. A 2-fold level was chosen as this corresponded to one standard deviation in a double 
determination of -1800 genes. Centromeres and heterochromatic regions were excluded from data analysis. 



grade I and II, respectively, tumors 733 and 827 were staged as pTI 
(invasive Into submucosa), 733 was staged as solid, and 827 was 
staged as papillary, both grade III. 

mRNA Prepa/at/on -Tissue biopsies, obtained fresh from surgery, 
were embedded Immediately In a sodlum-guanldlnlum thlocyanate 
solution and stored at -80 °C. Total RNA was Isolated using the 
RNAzoI B RNA isolation method (WAK-Chemle Medical GMBH). 
potiW* RNA was isolated by an oBgo(dT) selection step (Oligotex 
mRNA kit; Qiagen). 

cflrVA Preparation-! jxg of mRNA was used as starting material. 
The first and second strand cDNA synthesis was performed using the 
Superscript® choice system (invitrogen) according to the manufac- 
turer's Instructions but using an oligo(dT) primer containing a T7 RNA 
polymerase binding site. Labeled cRNA was prepared using the ME- 
QAscrlp® /V? vitro transcription kit (Amblon). Biotin-Iabeled CTP and 



UTP (Enzo) was used, together with unlabeled NTPs In the reaction. 
Following the In vitro transcription reaction, the unIrKX*porated nu- 
cleotides were removed using RNeasy columns (Qiagen). 

Array Hybridization and Scanning-Array hybridization and scan- 
ning was modified from a previous method (13). 10 & of cRNA was 
fragmented at 94 °C for 35 mln In buffer containing 40 mM Tris 
acetate, pH 8. 1 , 1 00 mM KOAc, 30 mM MgOAc. Prior to hybridization, 
the fragmented cRNA In a 6x SSPE-T hybridization buffer (1 m NaCI, 
10 mM Tris, pH 7.6, 0.005% Triton), was heated to 05 «C for 5 mln, 
subsequently cooled to 40 °C. and loaded onto the Affymelrix probe 
array cartridge. The probe array was then incubated for 18 h at 40 °C 
at constant rotation (BO rpm). The. probe array was exposed to 10 
washes In 6X SSPE-T at 25 °C followed by 4 washes in 0.5X SSPE-T 
at 50 °C. The blotlnylated cRNA was stained, with a streptavidin- 
phycoerythrin conjugate, 1 0 /ug/ml (Molecular Probes) In 8x SSPE-T 
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Fig. 1— -continued 



for 30 mto at 25 °C followed by 10 washes In 6x SSPE-Tat 25 The 
probe arrays were scanned at 560 nm using a confocal laser scanning 
microscope {made for Affymetrix by Hewlett-Packard). The readings 
from the quantitative scanning were analyzed by Affymetrix gene 
expression analysis software. 

Mfcnosatefftte Ana/ysfe—Microsateitite Analysis was performed as 
described previously (14). Microsatellttes were selected by use of 
www.nobD^;nlh.gbv/genemap98 v and primer sequences were ob- 
tained from the genome data base at www.gcfc.org. DMA was extracted 
from tumor and blood and amplified by PCR in a volume of 2 0 jd for 35 
. cycles. The amplfcons were denatured and electrophoresed (or 3 h In an 
ABI Prism 377. Data were collected In (he Gene Scan program for 
fragment analysis. Loss of heterozygosity was defined as less than 33% 
of one allele detected In tumor ampncons compared with blood. 

Proteomtc Anafysls— TCCs were minced Into small pieces and 
homogenized in a small glass homogenlzer In. 0.5 ml of lysis solution. 
Samples were stored at -20 °C until use. The procedure for 2D gel 
electrophoresis has been described In detail elsewhere (15, 16). Gels 
were stained with sliver nitrate and/or Coornassle Brilliant Blue. Pro- 
teins were Identified by a combination of procedures that Included 
micros equendng, mass spectrometry, two-dimensional gel Western 
Immunoblottlng. and comparison with the master two-dimensional gel 
Image of human kerat/nocyte proteins; see Wobase.dk/cgRJn/cens, 

CGW— Hybridization of differentially labeled tumor and normal DNA 
to normal metaphase chromosomes was performed as described 
previously (10). Ruoresceln-labeled tumor DNA (200 ng), Texas Red- 



labeled reference DMA (200 ng), and human Cot-1 DNA (20 pg) were 
denatured at 37 °C for 5 mln and applied to denatured normal met- 
aphase slides. Hybridization was at 37 °C for 2 days. After washing, 
trie slides were counters tain ed with 0.1 5 pg/ml 4>6-cfiamldino-2-p he- 
nytmdole In an anti-fade solution. A second hybridization was per- 
formed for all tumor samples using fluorescein- labeled reference DNA 
and Texas Red-labeled tumor DNA (Inverse labeling) to confirm the 
aberrations detected during the initial hybridization. Each CGH ex- 
periment also included a normal control hybridization using fluores- 
cein- and Texas Red-labeled normal DNA. Digital image analysis was 
used to Identify chromosomal regions with abnormal fluorescence 
ratios/ indicating regions of DNA gains and tosses. The average 
greentred fluorescence intensity ratio proffles were calculated using 
four Images of each chromosome (eight chromosomes total) with 
normalization of the greenrred fluorescence Intensity ratio for the 
entire metaphase and background correction. Chromosome identifi- 
cation was performed based on 4.6-dlamldIno-2-phenyUndole band- 
ing patterns. Only Images showing uniform high Intensity fluores- 
cence with minimal background staining were analyzed. All 
centromeres, p arms of acrocentric chromosomes, and heterochro- 
malic regions were excluded from the analysis. 

RESULTS 

Comparative Genomic Hybridization— The CGH analysis 
identified a number, of chromosomal game and losses In the 
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Table I 

Correlation between aJterotlohs detected by COM and by expression monitoring 

Top, CGH used as independent variable (If CGK alteration - what expression ratio was found); bottom, altered expression used. as 
Independent variable fcf expression alteration - what CGH deviation was found). 



CGH alterations 



Tumor 733 vs. 335 
Expression change clusters 



Concordance CGH alterations 



Tumor 827 vs. 532 
Expression change clusters 



Concordance 



13 Gain 10 Up-reguJaUcn 

0 Down-regulation 

3 No change 
10 loss 1 Up-reguiation 

S Powrwegulatlon 

4 No change 



77% 
60% 



10 Gain. 8 Up-regulation 

0 Oowiwegulatlon 
2 No change 

12 Loss 3 Up-regulation 

2 Down regulation 
7 No change 



Expression change dusters 



Tumor 733 vs. 335 
CGH alterations 



Concordance Expression change dusters 



Tumor 827 vs. 532 
CGH alterations 



8096 

17% 

Concordance 



16 Up-regulaiion 
21 Down-regulation 
15 No change 



11 Gain 

2 Loss 

3 No change 
1 Gain 

8 Loss 

12 No change 
3 Gain 

3 Loss 

9 No change 



69% 
38% 
60% 



17 Up-regulafion 
9 Dovvn^regulation 
21 No change 



10 Gain 

5 Loss 

2 No change 
OGaln 

3 Loss 

6 No change 
1 Gain 

3 Loss 

17 No change 



59% 
33% 
81% 



' two Invasive tumors (stajge pT1 , TCCs 733 and 827), whereas 
the two non-invasive papillomas (stage pTa, TCCs 335 and 
532) showed only 9p-, 9q22-q33-, and X-, and 7+, 9q-, 
and Y-, respectively- Both invasive tumors showec! changes 
(1q22-24-h, 2q14.1-qter-, 3q12-q13.3-, 6q12-q22-, 
9q34+, 11q12-q13+, 17+, and 20q11.2-q12+) that are typ- 
ical for their disease stage, as well as. additional alterations, 
some of which are shown in Fig. 1. Areas with gains and 
losses deviated from the normal copy number to some extent, 
and the^average numerical deviation from normal was 0.4-fold 
In the case of TCC 733 and 0.3-fold for TCC 827. The largest 
changes, amounting to at least a doubling of chromosomal 
content, were observed at 1q23 in TCC 733 (Fig. 1<4) and 
20q12 in TCC 827 (Fig. 18). 

mRNA Expression in Relation to DNA Copy Number—The 
mRN/yevets from the' two irweslvejtuniorjs, fTCCs 827 and 
733) were compared with the two non-invasive counterparts 
(TCCs 532 and 335). This was done In two separate experi- 
ments In which we compared TCCs 733 to 335. and 827 to 
532, respectively, using two different scaling settings for the 
arrays to rule out scaling as a confounding parameter. Ap- 
proximately 1 ,800 genes that yielded a signal on the arrays 
were searched In the Unigene and Genemap data bases for 
chromosomal location, and those with a known location 
(1096) were plotted as bars covering their purported locus. In 
that way it was possible to construct a graphic presentation of 
DNA copy number and relative mRNA levels along the Indi- 
vidual chromosomes (Fig. 1). 

For each mRNA a ratio was calculated between the level In 
the Invasive versus the non-Invasive counterpart Bare, which 
represent chromosomal location of a gene, were color-coded 
according to the expression ratio, and onty differences larger 



than 2-fold were regarded as Informative (Rg. 1). The density 
of genes along the chromosomes varied, and areas contain- 
ing only one gene were excluded from the calculations. The 
resolution of the CGH method is very low, and soma of the 
outlier data may be because of the fact that the boundaries of 
the chromosomal aberrations are not known at high resolution. 

Two sets of calculations were made from the data For the 
first set we used CGH alterations as the independent variable 
and estimated the frequency of expression alterations in these 
chromosomal areas. In general, areas with a strong gain of 
chromosomal material contained a cluster of genes having 
Increased mRNA expression. For example, both chromo- 
somes 1q21-q25, 2p and 9q, showed a relative gain of more 
than 100% in DNA copy number that was accompanied by 
increased mRNA expression levels In the two tumor pairs (Fig. 
. 1). In most cases, chromosomal gains detected by CGH were 
accompanied by an Increased level of transcripts In both 
TCCs 733 (77%) and 827 (80%) (Table I, top). Chromosomal 
losses, on the other hand, were not accompanied by de- 
' creased expression in several cases, and were often. regis- 
tered as having unaltered RNA levels (Table I, top). The inabil- 
ity to detect RNA expression changes In these cases was not 
because of fewer genes mapping to the lost regions (data not 
shown). 

In the second set of calculations we selected expression 
alterations above 2-fold as the independent variable and es- 
timated the frequency of CGH alterations In these areas. As 
above, we found that Increased transcript expression corre- 
lated with gain of chromosomal material (TCC 733, 69% and 
TCC 827, 59%), whereas reduced expression was often de- 
tected In areas with unaltered CGH ratios (Table S, bottom). 
Furthermore, as a control we looked at areas with no atter- 
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3,6 



Mi* 



41 
AAA 
A 



' A 
AAA 



AAA, 
AA 



Expression changes 
' detected 



Expression changes 



Expression changes 
detected 



Expression changes 
not detected 



Tumor 827 versus 532 ' Tumor 733 versus 335 

Fto. 2. Correlation between maximum OGH aberration and the eblfity to detect expression change by oligonucleotide array 
monitoring. The aberration is shown as a numerical -fold change In ratio between Invasive tumors 827 (A) and 733 (4 ) arid their noninvasive 
counterparts 532 and 335. The expression change was taken from the Expression tine to the right In Fig. 1, which depicts the resulting 
expression change for a given chromosomal region. At least half of the mRNAs from a given region have to be either up- or down-regulated 
to be scored as an expression change. Ail chromosomal arms In which the CGH ratio plus or minus one standard deviation was outside the' 
ratio value of one were Included. ' 



atlon in expression. No alteration was detected. by CGH In. 
most of these areas (TCC 733, 60% and TOO 827, 81%; see 
Table I, bottom). Because the ability to observe reduced or 
increased mRNA expression clustering to a certain chromo- 
somal area dearly reflected the extent of copy number 
changes, we plotted the maximum CGH aberrations In the 
regions showing CGH changes against the ability to detect a 
change in mRNA expression as monitored by the oligonucleo- 
tide arrays (Fig. 2)v For both tumors TCC 733 (p < 0.015) and 
TCC 827 <p < 0.00003) a highly significant correlation was 
observed between the level of CGH ratio change (reflecting 
the DNA copy number) and alteratTons detected by the array 
based technology "(Rg. 2fc Similar data were obtained when 
areas with altered expression were used as Independent vari- 
ables. These 'areas correlated best with CGH when the CGH 
ratio deviated 1 .6- to 2.0-fold (Table I. bottom) but mostly did 
not at lower CGH deviations. These data probably reflect that 
loss of an allele may only lead to a 50% reduction In expres- 
sion level, which Is at the cut-off point for detection of expres- 
sion alterations. Gain of chromosomal material can occur to a 
much larger extent, 

Mlcrosatetllte-based Detection of Minor Areas of Loss- 
es—In TCC 733, several chromosomal areas exhibiting DNA 
amplification were preceded or followed by areas with a nor- 
mal CGH but reduced mRNA expression (see Fig. 1 , TCC 733 
chromosome 1q32, 2p21, and 7q21 and q32, 9q34, and 
10q22). To determine whether these results were because of 
undetected toss of chromosomal material in these regions or 



because of other non-structural mechanisms regulating tran- 
scription, we examined two mlcrosateilltes positioned at chro- 
mosome 1q26-32 and two at chromosome 2p22. Loss of 
heterozygosity (LOH) was found at both 1q25 and at 2p22 
Indicating that minor deleted areas were not detected with the 
resolution of CGH (Fig. 3). Additionally, chromosome 2p In 
TCC 733 showed a CGH pattern of gain/no change/gain of 
DNA that correlated with transcript increase/decrease/in- 
crease. Thus, for the areas showing increased expression 
there was a correlation with the DNA copy number alterations 
(Fig! 1 A). As indicated above, the mRNA decrease observed In 
the middle of the chromosomal -gain was because of LOH, 
implying that one of the mechanisms" for mRNA down-regu- 
lation may be regions that have undergone smaller tosses of 
chromosomal material. However, this cannot be detected with 
the resolution of the CGH method. 

In both TCC 733 and TCC 827, the telomeric end of chro- 
mosome 11p showed a normal ratio In the CGH analysis; 
however, clusters of five and three genes, respectively, lost 
their expression. Two microsateUHes (D11S1760, D11S922) 
positioned close to MUG2, IGF2, and cathepsln D Indicated 
LOH as the most likely mechanism behind the loss, of expres- 
sion (data not shown). 

. A reduced expression of mRNA observed in TCC 733 at 
chromosomes 3q24» 11p11, 12p12.2, 12q21.1, and 16q24 
and In TCC 827 at chromosome 11p15.5, 12p11, 15q11.2, 
and 18q12 was also examined for chromosomal losses using 
mlcrosateilltes positioned as close as possible to the gene loci 
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Fta. 3. Mlcrosatellite analysis of loss of heterozygosity. Tumor 
733 showing loss of heterozygosity at chromosome 1q25, detected 
(a) by D1S21S close to Hu class I histocompatibility antigen (gene 
number 3d .in Rg. 1), (b) by D1S2735 close to cathepsin E (gene 
number 41 In Fig. U and (c) at chromosome 2p23 byJJ2§^25^close 
to general 0-spectrin (gene number 1 1 on Fig. 1) and of W) tumor 827 
showing toss at heterozygosity at chromosome 18q12 by S18S1118 
dose to mitochondrial 3-oxoacyl-coenzyme A thiol ase (gene number 
12 In Fig. 1). The upper curves show the electropherogram obtained 
from normal DMA from leukocytes (W), and the tower curves show the 
electropherogram from tumor DMA (7). In all cases one allele Is 
partially lost In the tumor ampllcon. 

showing reduced mRNA transcripts. Only the mlcrosatelilte 
positioned at 18q12 showed LOH (Rg. 3), suggesting that 
transcriptlonal.down-regulatlon of genes In the other regions 
may be controlled by other mechanisms. 

Relation between Changes In mRNA and Protein Levels- 
2D-PAGE analysis, In combination with Coomassle Brilliant 
Blue and/or silver staining, was carried out on all four turnors 
using fresh biopsy material. 40 wail resolved abundant known 
proteins migrating In areas away from the edges of the pH 



0.1 



B 

A 



F 

• t i» 



•jc 



A B 



Reduced 
protein 



Unaltered 
protein 



increased 
protein 



Rg. 4. Correlation between protein levels as Judged by 20* 
PAGE and transcript ratio. For comparison proteins were divided In 
three groups, unaltered In level or up* or down-regulated IjwrizontaJ 
axis). The mRNA ratio as determined by oligonucleotide arrays was 
plotted for each gene {vertical axis), a, mRNAs that were scored as 
present In both tumors used for the ratio calculation; A, mRNAs that 
were scored as absent in the Invasive tumors (along horizontal axis) or 
as absent in* noninvasive reference {top of figure). Two different 
scafings were used to exclude scafing as a con founder, TCCs 827 
and 532 (AA) were scaled with background suppression, and TCCs 
733 and 335 (#C) were scaled without suppression. Both compari- 
sons showed highly significant ip < 0.005) differences' In mRNA ratios 
between the groups. Proteins shown were as follows: Group A (from 
loft), phosphoglucomutase 1 , glutathione transferase class p number 
4, fatty acid-binding protein homotogue, cytokeratin 15, and cyto- 
keratin 13; B (from /etf), fatty; acld-Wnding protein homotogue, 28-kDa 
heat shock protein, cytokeratin 13, and calcycOn; C<from fen), a-eno- 
lase, hnRNP B1, 28-kDa heat shock protein, 14-3-3-c and 
pre-mRNA splicing factor; D t mesothelfal keratin K7 (type II); B (from 
top), glutathione S -transferase and mesotheOal keratin K7 (type. II); 
F (from top and tefi), adenylyl cyclase-assoclated protein, E-cadherin, 
keratin 19, calglzzarin, pnosphogJycerate mutase, annaxJn IV, cy- 
toskeletaJ y-actln, hnRNP At, Integral membrane protein calnexfn 
flP90), hnRNP H, brain-type cfathrin Ught chaln-a, hnRlsP F, 70-kDa 
heat shock protein, heterogeneous nuclear ribonudeoproteJn A/B, 
translattonally controlled tumor protein, liver gtyceraktehyde-3-phos- 
phate dehydrogenase, keratin 8, aldehyde reductase, and Na,K- 
ATPase 0-1 sub unit; G, (from fop and tefl}, TCP20, calglzzarln, 70- 
kDa heat shock protein, calnexln, hnRNP H, cytokeratin 15, ATP' 
synthase, keratin 19,.trlosephosphate Isomerase, hnRNP F f liver give- 
eraldehyda-3-phosphatase dehydrogenase, glutathione S-tranafer- 
ase-ir, and keratin 8; H (from fen), plasma gelsolin, autoantigen cai- 
retlcuBn, thloredoxin, and NAD+ -dependent 15 hyoVoxyprostaglandln 
dehydrogenase: / (from top), prolyl 4-hydroxylase 0-subuntt, cyto- 
keratin 20, cytokeratin 17, prohibition, and fructose 1,6-blphos*- 
phatase; J annexln II; K, ennexin IV; L (from top and fefl), 90-kDa heat 
shock protein, prolyl 4 -hydroxylase 0-subunlt, a-enolase, GRP 78, 
cyciophilln, and cofilln. 

gradient, and having a known chromosomal location, were 
selected for analysis In the TCC pair 827/532. Proteins were 
Identified by a combination of methods (see "Experimental 
Procedures"). In general there was a highly significant corre- 
lation ip < 0.005) between mRNA and protein alterations (Fig. 
4). Only one gene showed disagreement between transcript 
alteration and protein alteration. Except for a group of cyto- 
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keratins encoded by genes on chromosome 17 (Fig. 5) the 
analyzed proteins did not belong to a particular family. 28 well 
focused proteins whose genes had a know chromosomal 
location were detected in TCCs 733 and 335, and of these 19 
correlated (p < 0.005) with the mRNA changes detected using 
the arrays (Fig. 4). For example, PA-FABP was highly ex- 
pressed In the non-invasive TCC 335 but lost In the Invasive 
counterpart (TCC 733; see Rg. 5). The smaller number of 
proteins detected In both 733 and 335 was because of the 
smaller size of the biopsies that were available. 

11 chromosomal regions where CGH showed aberrations 
that corresponded to the changes in transcript levels also 
showed corresponding changes In the protein level (Table ||). 
These regions Included genes that encode proteins that are 
found to be frequently altered In bladder cancer, namely 
cytokeratins 17 and 20, annexins II and IV, and the fatty 
acid-binding proteins PA-FABP and FBP1. Four of these pro- 
teins were encoded by genes In chromosome 17q, a fre- 
quently amplified chromosomal area in invasive bladder 
cancers. 

DISCUSSION 

Most human cancers have abnormal DNA content, having 
lost some chromosornal parts and gained others. The present 
study provides some evidence as to the effect of these gains 
and losses on gene expression in two pairs of noninvasive 
and Invasive TCCs using high throughput expression arrays 
and proteomJcs, In combination with COM In general, the 
results showed that there ts a clear individual regulation of the 
mRNA expression of single genes, which In some cases was 
superimposed by a DNA copy number effect Inmost cases, 
genes located In chromosomal areas with gains often exhib- 
ited increased mRNA expression, whereas areas showing 
losses showed either no change or a reduced fnRNA expres- 
sion] The latter might be because of the fact that losses most 
often are restricted to loss of one allele, and the cut-off point ' 
for detection of expression alterations was a 2-fold change, 
thus being at the border of .detection. In several cases, how- 

.Table II 



Proteins whose expression level correlates with both mRNA and gene dose changes 


Protein 


Chromosomal location 


Tumor TGC 


CGH alteration 


Transcript alteration" 


Protein alteration 


Anrtexln II 


1q21 


733 


Gain 


Abs to Pres" 


Increase 


Annexln IV 


2p13 


733 


Gain 


3.9-Fold up 


Increase 


Cytokeratin 17 


17q12-q21 


827 


Gain . 


3.8-Fold up 


increase 


Cytokeratin 20 


17q21.1 


827 


Gain 


53-Fold up 


Increase 


(PA-)FABP 


8q21.2 


827 


Loss 


10-Fold down 


Decrease 


FBP1 


9q22 


627 


Gain 


2.3-Fold Up 


Increase 


Plasma getsofln 


9q31 


827 


Gain 


Abs to Pres 


Increase 


Heat shock protein 28 


15q12-q13 


827 


less 


23-Fold up 


Decrease 


Prohibffln 


17q21 


827/733 


Gain 


3.7-/2.5-Fold up* 


Increase 


Proryl-4-hydroxyl 


17q25 


827/733 


Gain * 


6.7-/1 ,6-Fotd up 


Increase 


hnRNPBI 


7p15 


827 


Loss 


2.5-FokJ down 


. Decrease 



* Abs, absent; Pres, present. 

» In cases .where the corresponding alterations were found in both TCCs 827 and 733 these are shown as 827/733, 




Fig. 5. Comparison of protein and transcript levels In Invasive 
and non-Invasive TCCs. The upper part of the figure shows a 2D gel 
(tofl) and the oligonucleotide array {righQ of TCC 532. The red rectan- 
gles on the upper get highlight the areas that are compared below. 
Identical areas of 2D gels of TCCs S32 and 827 are shown below. 
Clearly, cytokeratins 13 and 15 are strongly down-regulated In TCC 
827 (red annotation). The tile on the army containing probes for 
cytokeratin. 15 Is enlarged below the array (red arrow) from TCC 532 
and Is compared with TCC 827. The upper row of squaws in each tile 
corresponds to perfect match probes; the lower row corresponds to 
mismatch probes containing a mutation (used for correction for un- 
speciflc binding). Absence of signal ts depleted as black, and the 
higher the signal the lighter the cotor. A high transcript level was 
detected fn TCC 532 (6151 units) whereas a much tower level was 
detected In TCC 827 (absence of signals}. For cytokeratin 13, a high 
transcript - level was also present In TCC 532 (15659 units), and a 
much lower level was present in TCC 827 (623 units). The 2D gels at 
the bottom of the figure (left) show levels of PA-FABP and adipocyte- 
FABP in TCCs 335 and 733 (invasive), respectively. Both proteins are 
down-regulated In the invasive tumor. To the right we show the array 
tiles for the PA-FABP transcript. A medium transcript level was de- 
tected in the case of TCC 335 (1277 units) whereas very low levels 
were detected in TCC 733 (166 units). IEF, Isoelectric focjusing. . . , 
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ever, an Increase or decrease in DNA cop/ number was arm and that the use of cDNAmlcroarrays for analysis of DNA 
associated whh do novo occurrence or complete loss of tran- copy number changes will reach a resolution that cart resolve 
script, respectively. Some of these transcripts could not be these changes, as has recently been proposed. (2). The outlier 

. detected In the noninvasive tumor but were present si rela- data were not more frequent at the boundaries of the CGH 
tively high levels In areas with DNA amplifications irt the inva- aberrations. At present we do not know the mechanism be- 
slve tumors (e.g. In TCC 733 transcript fern cellular llgand of hind chromosomal aneuptoldy and cannot predict whether 
annexto 11 gene (chromosome 1q21) from absent to 2670 chromosomal gains will be transcribed to a larger extent than 
arbitrary Units; In TCC 827 transcript from small prollne-rlch the two native alleles. A mechanism as genetic Imprinting has 
protein 1 gene (chromosome 1q12-q21.l) from absent to an Impact on the expression level In normal cells and Is often 
i326 arbitrary units). It may be anticipated from these data reduced In tumors. However, the relation between Imprinting 
that significant clustering of genes with an Increased express and gain of chromosomal malerialis not known, 
slon to a certain chromosomal area indicates an Increased We regard It as astrength of this Investigation that we were 
likelihood of gain of chromosomal material In this area, able to compare Invasive tumors to benign tumors rather than 

Considering the many possible regulatory mechanisms act- to norjnal urotheliiim, as the tumors studied were biologically 
Ing at the level of transcription, It seems striking that the gene very close, and probably may represent successive steps In 
dose effects were so clearly detectable In gained areas. One the progression of bladder cancer. Despite me limited amount 
hypothetical explanation may Re In the loss of controlled of fresh tissue available ft was possible to a 
rnethylation In tumor cells (17-19). Thus, It may be possible . state of the art methods. The observed correlation between 
that in chromosomes with increased DNA copy numbers two DNA copy number and mRNA expression Is remarkable when 

. or more alleles could be demethylated simultaneously leading one considers that different pieces of the tumor biopsies were 
to a" higher transcription level, whereas In chromosomes with used for the different sets of experiments. This Indicate that 
tosses the remaining allele could be partly methylated, turning bladder tumors are relatively homogenous, a notion recently 
off the process (20) 21). A recent report has documented a supported by CGH and LOH data that showed a remarkable 
ptoldy regulation of gene expression In yeast, but In this case all similarity even between tumors and distant metastasis (10, 23). 
the genes were present In the same ratio £2), a situation that is In the few cases analyzed, .mRNA. and protein levels 
not analogous to the* of cancer cells, wlifch show marked showed a striking correspondence although In some cases 
crvbmosomal aberrations, as well as gene dosage effects.* we found discrepancies that may be attributed to trartslational 
Several CGH studies of bladder cancer have shown that regulation, post-transiatlonal processing, protein degrada- 
some chromosomal aberrations ere common at certain tion, or a combination of these. Some transcripts belong to 
stages of disease progression, often occurring In more than 1 undertranslated mRNA pools, which are associated with few 
of3tumonUnpTatumors,.th^ translationally Inactive rlbosomes; these pools, however, 

(2. 6), and.ln;pT1 tumors. 2q-,11p- 1 1q-, 1q+, 5p+, 8q+, seem to be rare (24). Protein degradation, for example, may 
17q+, and 20q+ (2-4, 6, 7). The pTa tumors studied here be very Important. in the case of polypeptides with a short 
showed similar aberrations such as 9p- and 8q22-q33- and half-life (e-g. signaling proteins), A poor correlation between 
9q- and respectively. Ukewlse, the two minimal invasive mRNA and protein levels was found' in liver cells as deter- 
pT1 tumors showed aberrations that are commonly seen at mined by arrays and 2D-PAGE (25), and a moderate correla- 
that stage, and TCC 827 had a remarkable resemblance to the tior* was recently reported by Ideker etaL (26) In yeast 
commonly seen pattern ofOosses and gains, such as 1q2?-24 . Interestingly, our study revealed a much better correlation 
amplification (seen in both tumors}, 11q14-q22 loss, the latter between gained chromosomal areas and Increased mRNA 

. often linked to 17 q+ (both tumors), and 1q+ and 9p~, often levels than between loss of chromosomal areas and reduced 
linked to 20q+ and 11 q13+ (both tumors) (7-9). These ob- mRNA levels. In general, the level of CGH change determined 
servations indicate that the pairs of tumors used In this study the ability to detect a change in transcript, One possible 
exhibit chromosomal changes observed in many tumors, and explanation could be that by losing one allele the change In 
therefore the findings could be of general importance for mRNA level is not so dramatic as compared with gain of 
bladder cancer. material, which can be rather unlimited and may lead to a 

Considering that the mapping resolution of CGH Is of about severalfold increase In.gene copy number resulting In a much 
20 megabases it is only possible to get a crude picture of higher impact on transcript level, the latter would be much 
chromosomal Instability using this technique. Occasionally, easier to detect on the expression anrays as the cut-off point 
we observed reduced transcript levels close to or Inside re- was placed at a 2-fbld level so as not to be biased by noise on 
glons with increased copy numbers. Analysis of these regions the array. Construction of arrays with a better signal to noise 
by positioning heterozygous mlcrosatelRtes as close as pos« ratio may In the future allow detection of lesser than 2-fold 
slble to the locus showing reduced gene expression revealed alterations in transcript levels, a feature that may facilitate the 

• loss of heterozygosity In several cases. It seems likely that analysis of the effect of loss of chromosomal areas on tran- 
multiple and different events occur along each chromosomal script levels. 



44 Molecular & Cellular Profeomics 1. 1 



Gene Copy Numbers, Transcripts, and Protein Levels 



> In deven cases we found a significant correlation between 
DNA copy number, mRNA expression; and protein level. Four 
of these proteins were encoded by genes located at a fre- 
quently amplified area In chromosome 17q.. Whether DNA 
copy number Is one of the mechanisms behind alteration of 
these eleven protelns'is at present unknown and will have to 
be proved by other methods using a larger number of sam- 
ples. One factor making such studies complicated Is the large 
extent of protein modification that occurs after translation, 
requiring &nmunoldentlflcatlon and/or mass spectrometry to 
correctly Identify the proteins In the gels: 

In conclusion, the results presented In this study exemplify 
the large body of knowledge that may be possible to gather In 
the future by combining state of the art techniques that follow 
trie pathway from DNA to protein (26). Here, we used a. tradi- 
tional chromosomal CGH method, but in the future high reso- 
lution CGH based on mlctoarrays with many thousand radiation 
hybrid-mapped genes will Increase the resolution and informa- 
tion defivedfrom these types of experiments (2). Combined with 
expression arrays analyzing transcripts derived from genes with 
known locations, and 2D gel analysis to obtain Information at 
the post-transtatlonai level, a dearer and more developed un- 
derstanding of the tumor genome will be forthcoming. 
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